RIKEN-TH-75 
hep-th/0608031 

Dynamical Generation of Non-Abelian Gauge Group 
via the Improved Perturbation Theory 

T. Aoyama, T. Kuroki, and Y. Shibusa 
Theoretical Physics Laboratory, RIKEN, Wako, 351-0198, Japan 

Abstract 

It was suggested that the massive Yang-Mills-Chern-Simons matrix model has three phases and 
that in one of them a non-Abelian gauge symmetry is dynamically generated. The analysis was at 
the one-loop level around a classical solution of fuzzy sphere type. We obtain evidences that three 
phases are indeed realized as nonperturbative vacua by using the improved perturbation theory. It 
also gives a good example that even if we start from a trivial vacuum, the improved perturbation 
theory around it enables us to observe nontrivial vacua. 
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I. INTRODUCTION 



It is now widely believed that some gauge theories or matrix models would provide non- 
perturbative formulations of string theories in the large- N limit l[ [2], ^] . In such approaches, 
it is quite important to clarify how they quantize gravity, namely, how they reconcile the 
general relativity and quantum theory. It is also of great importance to confirm whether at a 
low energy regime they can be effectively described by quantum field theories of gauge fields 
with the non-Abelian gauge group of the standard model type. However, in the context 
of the matrix models it is not yet known how the gauge group of the standard model is 
dynamically generated. It is natural that it is embedded in the original U (N) gauge group 
of a matrix model we start from and that it emerges dynamically in the large- N limit 4j 
via the Higgs mechanism or the confinement, for example. Therefore, it is interesting and 
necessary to construct matrix models that incorporate such a mechanism and to study them 
in a nonperturbative manner. 

As a first step toward this program, the bosonic massive Yang-Mills-Chern-Simons 
YMCS) matrix model was analyzed in Ref. which is obtained by the large- N reduction 
6| of the U (N) Yang-Mills theory with the Chern-Simons term in three dimensions and by 
adding the mass term. Without the mass term, it was shown in Ref. {7] that the YMCS 
matrix model has two phases, namely, the Yang-Mills phase and the fuzzy sphere phase. 
In the former phase observables show qualitatively the same behavior as those in the pure 
Yang-Mills model without the Chern-Simons term, while in the latter a single fuzzy sphere 
is a stable and dominant configuration due to the Chern-Simons term. It was claimed in 
Ref. J5] that if the mass term is present there appears an interesting third phase ("5-brane 
phase") in addition to these two phases, in which a dominant configuration is expected to 
be O(N) coincident fuzzy spheres of small size and thus a non-Abelian gauge group with 
rank of O(N) is dynamically generated as a subgroup of the original U(N) gauge group 
Although in the fuzzy sphere phase comparison between the one-loop calculation around 
a single fuzzy sphere configuration and the Monte Carlo data shows one-loop dominance, 



namely higher 
are concerned 



oop corrections are suppressed in the large- N limit as long as 1PI diagrams 
, it is not a priori trivial whether the analysis at the one-loop level is 
sufficient even in the presence of mass term. 

In order to examine whether the 5-brane phase found in Ref. [5] is indeed realized as a 
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nonperturbative vacuum, we have to investigate the massive YMCS matrix model in a non- 
perturbative method. Toward this aim, we employ here the improved perturbation theory 
developed in Ref. [9], which is motivated by the improved mean field approximation (IMFA). 
This concept is also called the Gaussian approximation or Gaussian expansion method and 
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111 ], which has been further inves 



was introduced to a large- N reduced model in Ref. 

tigated in various _ flflQ. This technic i 8 a so rt of variations! .nethod |li 
to extract nonperturbative information on the vacuum of a theory. The IMFA is indeed 
reduced to the standard mean field approximation at the leading order, and the inclusion of 
higher-order contributions may be viewed as systematic improvement of the approximation. 



The IMFA has 
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y to a number of models to yield nonperturbative 
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201 ] . In particular, the IMFA was applied to the 



IIB matrix model first in Ref. [12j and provided the interesting nonperturbative observation 
that our four-dimensional space-time is realized in the IIB matrix model as the anisotropy of 
eigenvalue distribution of ten bosonic Hermitian matrices. This result was further confirmed 
in Refs. flflQQ- 

The improved perturbation theory was developed as such a reformulation of the IMFA 
that it is considered as a reorganization of a perturbative series. It opened a way to taking 
account of more general types of action than that had been treated in the IMFA or the 
Gaussian expansion method. Thus in the present article we examine the matrix model with 
Chern-Simons term that has a cubic interaction term in the framework of the improved 
perturbation theory. We expect that the observations obtained in this work would be in- 
structive to future studies of more nontrivial and interesting nonperturbative phenomena via 
the improved perturbation theory, in particular, dynamical generation of gauge symmetry 
in the IIB matrix model. 

In the prescription of the IMFA and the improved perturbation theory, a set of artificial 
parameters are introduced through a nominal deformation in the original action of the 
model. It leads to a sort of consistency conditions for these parameters that the original 
theory should not depend on them. In the space of these parameters, the conditions may 
be realized as a region in which physical quantities are least sensitive to the variation of the 
parameters [lj3]. Such a region is called "plateau". As stressed in Refs. jjj, formation 
of a plateau provides an essential criterion to see whether or not the improved perturbation 
theory works. It is shown in Ref. {20!] that in a system which exhibits a phase transition, 



3 



there may appear several plateaux. Therefore, it is interesting to discuss the phase structure 
of the massive YMCS matrix model in the framework of the improved perturbation theory 
through the patterns of emergence of plateaux that correspond to various phases. 

Another motivation of applying the improved perturbation theory to the massive YMCS 
matrix model comes from the fact that even if we start from the perturbative vacuum, 
the improved perturbation theory can provide information on a nonperturbative vacuum 
which does not always correspond to a classical configuration. Therefore, in principle, even 
if we start from the trivial vacuum which corresponds to the Yang-Mills phase, we would 
be able to identify the fuzzy sphere phase and the 5-brane phase when they are realized 
as nonperturbative vacua. We will see that this is indeed the case. This is in contrast 



with previous ana 
around them 5, 7 



yses 
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of fuzzy sphere type configurations based on the perturbation theory 



221. 



The organization of this paper is as follows: in the next section we review the massive 
YMCS matrix model and its three phases. In Section IHIl we also give a review of the 
improved perturbation theory, which is applied to the massive YMCS matrix model in 
Section HVl The results are presented in Section [V] Section |VT] is devoted to discussions. 
The tables are put collectively at the end of the paper. 

II. THE MODEL 

The massive Yang-Mills-Chern-Simons matrix model is defined by the action 

S = Nti f-~[Ai> M 2 + l -ae^ vX A^A u A x + \^ 2 A 2 }j , (2.1) 

where A^ (// = 1 ~ 3) are N x N Hermitian matrices. We will consider this model in the 
large- N limit. We can obtain this model by adding the mass term to the large- N reduced 
model of the U (N) Yang-Mills theory with the Chern-Simon term in three dimensions. This 
model has two parameters a, m and we are interested in the free energy of this model 

F = ~logZ, Z= I dA,e- s , (2.2) 



N 2 

as a function of them, by which we can determine the phase structure of this model. 
The equation of motion derived from Eq. (12. ip is 

[A x , [A\, Ap]] + % -ae^x[A v , A x ] + m 2 A^ = 0. (2.3) 



Apart from a trivial solution = 0, this model admits a classical solution of fuzzy sphere 
type given by 

A^ = X L^ n = 1 ~ 3, (2.4) 
where L^s is the iV-dimensional representation of SU{2) algebra satisfying 



(2.5) 



and x = ( a + Vet 2 — 8m 2 ) /A. If m 2 > a 2 /8, this type of solution ceases to exist. Gener- 
ically L M can be a reducible representation of the SU (2) algebra which consists of the rii- 
dimensional irreducible representations for i — 1 ~ s, where Yli=i n i = N. In this case 
without loss of generality A^ can be brought into the form 



/ r(ni) 



A, 



X 



\ 



\ 



(2.6) 



where denotes the n-dimensional irreducible representation. Substituting this into the 
action (I2.ip . we obtain the classical action of this solution as 



S* = Nf( X ) ~ 



m{n 2 - 1) 



i=l 



4 



(2.7) 



where f(x) = X 4 /2 — «X 3 /3 + m 2 % 2 /2. When a 2 > 9m 2 , /(%) is always negative, thus 
Eq. (12.71) implies that a single fuzzy sphere (s — 1, n\ — N) is the most dominant config- 
uration as long as a 2 and m 2 are not so small that S c \ ~ 0(1). Then the system is in the 
fuzzy sphere phase. On the other hand, when a 2 < 9m 2 , f(x) > and the trivial solution 
An = becomes stable. Moreover, whenever a is so small, the one-loop analysis around a 



background 



Af, = diag(x 



(1) T (2) 



5 , 



(2.8) 



shows that in general there is an attractive force between the eigenvalues of A„, and it makes 



their distribution shrink until the perturbative calculation becomes no longer valid 



23] 



Then they form a "solid ball" and the system is in the Yang-Mills phase. However, as shown 



in Ref. 



5j, this is not the end of the story. So far we have discussed the phase of the model 



based on the classical action (12 .7p . If we take account of the one-loop contribution to the free 



energy, we can find the third phase for the moderate a and m satisfying 8m 2 < a 2 < 9m 2 
fa. which OW copies of the smaU fuzzy sphere appear as the true vacuum. This phase 
is called "5-brane phase" in Ref. [5j because this configuration is analogous to the on e 
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that is interpreted [24[ as coinciding transverse 5-branes in the BMN matrix model 
By considering configurations with fuzzy spheres of various sizes, the phase diagram was 
proposed in Ref. \^ at the one-loop level. 1 The one- loop calculation is reliable in some 
parameter regions provided that the perturbative series around classical configurations are 
well-defined. 2 

However, in a general parameter region it is not clear that the one-loop calculation is 
sufficient and even if so, the perturbative series would be asymptotic. Therefore, it is 
interesting to confirm that the phase diagram in Ref. [5| is true even at a nonperturbative 
level by a totally different method. As such, we employ the improved perturbation theory 
developed in Ref. which works well even for some asymptotic series. 



III. THE IMPROVED PERTURBATION THEORY 

In order to describe the idea of the improved perturbation theory, we consider as an 
illustration a one-matrix model defined by the action 

S=ftr0 4 , (3.1) 

where is an N x N Hermitian matrix. Suppose we are interested in the free energy 

N 2 F = - log j d0exp ^-^tr0 4 ^ , (3.2) 

which cannot be computed by the standard perturbation theory for lack of the mass term 
(quadratic term) in Eq. (13. ip . Therefore, let us formally add and subtract the mass term in 
the action, one of which is regarded as a perturbation. Then we introduce a formal coupling 
constant g and construct a perturbative series in terms of g. g will be taken to be 1 in the 
end: 



2 2 



N F = - log / rf^exp ( -Mr ( |0 4 - g^4> 2 + , 



(3.3) 

9=1 



1 a and p in Ref. [5(] are written as a/2 and 2m /a respectively bousing the parameters in Eq. (|2.ip . 

2 In particular, the presence of the 5-brane phase is shown in Ref. [5j in the region a»l where the one-loop 
calculation is justified. 
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where mo is an artificially introduced parameter. It can be restated in the following way. If 
we define the free energy of a massive theory as 

N 2 F{m 2 ) = - log J d<\) exp ^-JVtr + ^0 2 ^ ^ , (3.4) 

the prescription above is equivalent to calculating the free energy of the massless theory 
(13. 2p from the massive theory with mass parameter replaced as 

^(0)1^=^0-^0)1,=!- (3-5) 

The perturbative series in terms of g obtained from the RHS of Eq. (13.51) is called the 
improved perturbative series. Note that if we can sum up this series at the full order, the 
free energy is completely independent of mo which is introduced artificially in order to make 
the perturbation theory feasible. 

In practice, we can calculate the RHS in Eq. (13.51) only up to a finite order. Then the 
dependence on mo emerges in the improved perturbative series. In order to obtain the exact 
value of the free energy, we need to determine the value of mo somehow. Here we adopt 



the principle of minimal sensitivity 15] as a guiding principle: the improved series should 
depend least on m . It is because the original theory does not depend on m which was 
introduced through a nominal shift of parameter. If there exists a region in the parameter 
space of m where the dependence on m vanishes effectively, the exact value would be 
reproduced there. We call such a region as "plateau". The principle of minimal sensitivity 
is realized as the emergence of plateau, in which the improved series stays stable against 
any variation of the artificial parameter. 

As for the one- matrix model (13.21) . we can identify a clear plateau even up to the 8th 
order, and there the improved perturbative series reproduces the exact value of the free 
energy with more than 99.5% accuracy [l3|. It is a striking result that although the ori gina l 



perturbative series of the massive theory has a radius of convergence g/m 4 < 1/12 26], 
we have obtained a good approximate value at m = by this technique. The plausibility 
and applicability of the improved perturbation theory is not yet clear in the mathematically 
strict sense. Nevertheless it has been successfully used in various models. Some examples 
are also given in Ref. 9}. 

Note that at the first order in g if we require dF[ mp /dm = as a "plateau" condition, 



this reduces to the self-consistency condition in the standard mean field approximation. In 



this sense, the improved perturbation series is a natural generalization of the mean field 
approximation. 

The scheme presented so far can be generalized to a generic series that may have more 
than one parameter. Thus we arrive at the concept of the improved Taylor expansion 9| 
as follows. Let us assume that an observable of a theory would be exactly described by a 
function F(X, £). Here A is a coupling constant and £ collectively represents parameters of 
the model such as a mass. Perturbation theory provides an expansion of F as a power series 
of A about A = 0, with nth coefficient denoted as 

oo 

F(A,0 = 5^A B / B (0. (3.6) 

n=0 

Next we consider a modification of the series according to the following prescriptions. First 
we perform a shift of parameters: 

A — > g\, 

£^£o + <?(£-£o), (3.7) 

where we have introduced g as a formal expansion parameter, and £o as a set of artificial 
parameters. We deform the series by the substitution (13.71) . and then we reorganize the 
series in terms of g up to g k and drop the 0(g k+l ) terms, and finally set g to 1. Thus we 
obtain the improved perturbative series i^ mp as 

F(A,0 — Fr(\,Z;Zo) = F(g\,Z + g(£ - Z ))\ upto 0(fffc) , fl=1 . (3.8) 

Here, a notation | up to , fc . represents the operation that we disregard the 0{g k+l ) terms 
and then put g to 1. 

More concretely, the improved perturbative series up to kth order F™ v is made by the 
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procedure: 

k 

F k (£) = ^A"/n(0 (3.9) 

n=0 
k 

n=0 
oo min(fc,£) 

= 5> £ E A» 7J — n (e /i^ceo) (3.H) 

£=0 n=0 ^ 

£=0 n=0 ^ >' 
k k—n 

— ^ E^E/f^-^/^o) (3-13) 

n=0 £=0 

= Fr p (A,e;e ). (3.14) 

Here, fjp is the fth derivative of / n (£) with respect to £. At first we apply the transformation 
( 13. 7p to the original fcth order perturbative series in Eq. (I3.10p . Next the series is then 
reorganized in terms of g in Eq. (13. lip . At this stage, we must take into account the extra 
g-dependence as £o + — Co) i n the coefficient f n . Then we drop the 0(g k+1 ) terms in 
Eq. (I3.12p . and finally set g to 1 in Eq. (13.131) . In this way we obtain the improved series 
F k mp (Ell. 

It turns out that by the procedure of the improved perturbation theory, each coefficient 
of the original series /„(£) are replaced by a particular combination of the coefficients in its 
Taylor expansion about a shifted point £o, 

k—n 1 

fn(0 — fn&Zo) = E T\ (3.15) 

e=o *" 

From this expression, we immediately find that even if Fk(£) in Eq. (I3.9P shows worse 
behavior, the improved perturbative series -F^ mp (£; £o) in Eq. (13.14j) often becomes mild as a 
function of £. 

In this procedure, the foregoing example of one matrix model corresponds to a case of 
A = 1 and £ = m 2 = where the improved perturbative series F™ v boils down to the 
expression (13.51) . 

Because the emergence of a plateau is an essential criterion to see whether the improved 
perturbation theory works well, a main task is identification of plateau with respect to the 
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artificial parameter £o at each value of the parameter £. Although we do not yet have a 
rigorous definition of plateaux, we could identify them by their properties investigated in 
|9|, [13|, [19|, |20[] as follows. The ideal realization of the plateau may have such a property that 



the improved series is totally independent of the artificial parameters £o i n a certain region. 
It involves a situation that all orders of derivatives of the improved series with respect to 
£o are zero in that region. However, such an ideal plateau is not realized in practical cases 
because we have series only of finite order. Typical profile of the improved series that forms 
plateau exhibits a flat region in which the series fluctuates bit by bit; it would consist of a 
number of local maxima and minima. Thus, we consider the accumulations of extrema as 
indications of plateau (or its candidates). 

In some cases it occurs that there is a region in which i^ mp becomes stable but varies 
gently without forming extrema. We should also take into account this region as a plateau. 
However, we have to mention that we often exclude such an asymptotic behavior that the 
series becomes flat at large values of parameters £o- 

For each value of £, there may appear more than one plateau, each of which corresponds 
to either a stable physical state or an unstable or metastable state (local minimum). By 
comparing the values of the improved free energy at these plateaux in the parameter space 
of £oj we are a ble to discuss which one is realized as a vacuum at the specified values of 
parameters £. The phase transition may also be argued on this basis in such a way that the 
state corresponding to the vacuum changes under the variation of the parameters £. 

It should be noted that although the free energy is a generating function of all observables 
which are obtained as derivatives with respect to corresponding parameters £, analyses of the 
improved perturbation theory should be applied independently to each observable of interest. 
This is because the prescription of improvement does not commute with differentiation with 
respect to the parameters £ and thus the formation of plateau may occur in a different 
manner by observables especially at lower orders of perturbation. 



IV. APPLICATION TO THE MASSIVE YANG-MILLS CHERN-SIMONS MA- 
TRIX MODEL 

In this section we discuss how to apply the improved perturbation theory to the massive 
YMCS matrix model. In order to determine the phase structure, we concentrate on the free 
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energy F(a,m 2 ) given in Eq. (12.21) and the expectation value of the Chern-Simons term 

A(a,m 2 ) = dF{a d f ] = (±tr (^V^)). (4- 1 ) 

as functions of two parameters a, m appearing in Eq. (12. ip . The reason we consider A(a, m?) 
is that it would be an order parameter of our model. Namely, discussion in SectionHTlsuggests 
that both F(a, m 2 ) and A(a, m 2 ) would be of 0(1) or even smaller in the Yang-Mills phase, 
while they would become of 0(1), and of 0(N 2 ) in the 5-brane phase, and the fuzzy sphere 
phase respectively, corresponding to the fuzzy sphere type configuration realized in each 
phase. 

According to the general prescription given in the previous section, what we should do 
is first to introduce a coupling constant A, which is set to be 1 in the end, into the original 
action (12. ip in order to define a perturbative series as 



S = Nti A u } 2 + -VXae^xA^AvAx + \rn 2 A 2 ^j 



(4.2) 



then to calculate F(a, m 2 ) or A(a, m 2 ) up to a certain order of A. Since we are interested in 
the large- N limit of our model, we fix the 't Hooft coupling A and consider only the planar 
diagrams. 3 Next we make the replacement 

A — ► gX, 

a — > a + g(oi — «o), 
m 2 — > ml + g{m 2 — m 2 ,) (4.3) 

in F(a, m 2 ) and A(a, m 2 ) to construct the improved perturbative series F imp (a, m 2 ; ao, tHq) 
and A irnp (a, m 2 ; ao, ttlq) for them. By searching for the plateau with respect to both ao and 
m\ introduced artificially, we can determine the true values of F(a, m 2 ) and A(a, m 2 ) for 
each a and m 2 . 

In order to obtain the free energy, we have to count and calculate all diagrams up to 
an order we hope. In fact, there are many diagrams which contribute to the free energy 
even in lower orders. In order to avoid this difficulty, we follow the prescription given in 
Ref. js]. Namely, the ordinary free energy is given by the Legendre transformation of the 



3 In other words, we take the N — > oo limit in such a way that the couplingconstants of four-point vertex 
and Chern-Simons term scale like Eq. (|4. 2|) which is the same as in Refs. 0). 
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2PI free energy given by the sum of diagrams in which there are no self-energy part and 
hence all propagators are the exact propagators. 4 In general, the number of diagrams which 
contribute to the 2PI free energy is drastically reduced compared to that of the ordinary 
free energy. This fact enables us to calculate the free energy in higher orders. In fact, we 
have performed it to the fifth order with respect to A. 

For the purpose of calculating the 2PI free energy, we assume the form of the exact 
propagator as 

{A^A V \) = CS^Sj". (4.4) 
Then C is the conjugate variable to m 2 in the sense of the Legendre transformation: 

OF 3^ 8G 3 2 

= ^ dC = ~2 m ' (45) 

where G is the 2PI free energy which is a function of a and C: G = G(a, C), from which 
we can obtain the free energy F = F(a,m 2 ) by the Legendre transformation with respect 
to C. 

Here we make a remark on the assumption (I4.4p . If we allow the most generic form of 
the exact propagator 

< Vi Vi) = <V/,, (4.6) 

it has 6N 4 parameters which become infinite in the large- N limit. Correspondingly, we have 

infinite number of parameters conjugate to them in the free energy. This is equivalent to 

introducing the most generic quadratic term in the action as a mean field. Evidently this 

seems intractable. Therefore we follow the same approach as in Refs. Namely, we 

impose appropriate symmetries on the exact propagator and assume its most generic form 

that is invariant under them. Then the number of the parameters in the exact propagator can 

be reduced to tractable one, and the free energy is written as a function of parameters allowed 

by the symmetries we have assumed. Thus we can find the most dominant configuration 

that respects them. In the present case let us assume the SO (3) Lorentz symmetry and 

the U(n) x U(n) x ■ ■ • U(n) xS m symmetry with N = nm. The latter symmetry reflects a 

v f\ 

"clustering" configuration [4] in which there are m clusters composed of n eigenvalues and 

hence the U (n) gauge symmetry is expected to be generated. S m is the permutation of 



The relation between the 2PI free energy and the Schwinger-Dyson equation is clarified in Ref. 
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these m clusters corresponding to the diffeomorphism on the fuzzy sphere. More precisely, 
the action of an element of S m on is defined by 

A H"')W) ^HWWvW)) for ° e Sm, ( ^ 4 ' 7 ' ) 
where (W) denotes the index i of the i'-th U(n) group in U(n) x U(n) x • • ■ U(n). Then it 

V v ' 

m 

is easy to see that the exact propagator is restricted to the form 

C ^ {U,) m [kk '\iv) = ( C5 n» + 05^5^)5^5/6^,5/', not summed on i',f. (4.8) 

From this expression we see that contributions of diagrams containing the second term are 
in general suppressed by 1/m. Therefore as long as m is O(N), we do not have to take them 
into account in the large- iV limit and the exact propagator can be given by Eq. (I4.4p . In 
fact, this is the case in the fuzzy sphere phase in which there exists a single fuzzy sphere 
and hence n — 1, m = N. On the other hand, the propagator should take the form in 
Eq. (14.41) also in the Yang-Mills phase, because the original U (N) gauge symmetry remains 
intact. However, in the 5-brane phase it is expected that there exist O(N) coincident small 
fuzzy spheres, which corresponds to the case n ~ O(N), m ~ 0(1). Therefore, from the 
viewpoint in Refs. we have to include the parameter D in Eq. (14.81) in calculating the 

2PI free energy, then the free energy would become a function of 3 parameters including one 
conjugate to D other than a and m. In order to avoid such a situation, we simply assume 
that the contribution of D in Eq. (14.81) would be irrelevant and that even in the 5-brane 
phase the exact propagator would take the form in Eq. (14. 4ft . Then we check consistency 
that the free energy F(a, m 2 ) and the expectation value of the Chern-Simons term A(a, m?) 
derived from G(a, C) show behavior expected from the configuration in the 5-brane phase. 
Another reason we anticipate (I4.4p even in the 5-brane phase is that from the standpoint 
of the noncommutative field theory, the original gauge symmetry is not broken, but is 
reinterpreted as the noncommutative gauge symmetry as shown in Refs. {§], [^j. In this 
sense, it is reasonable to assume (I4.4p in the Yang-Mills, 5-brane, and fuzzy sphere phases. 

Moreover, it is worth noticing that the improved perturbation theory with the ansatz in 
(I4.4p can investigate various vacua which are not restricted to three classical vacua as above 
as long as they satisfy the ansatz. 5 



5 However, the ansatz in (|4.4[) does not take account of configurations considered in Ref. Q where fuzzy 
spheres of various sizes appear. 
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and its 2-d plot in 



FIG. 1: Distribution of extrema of F imp (a, m 2 ; olq, itiq) with respect to tuq and ao at a = 0.0 
and m 2 = 0.0. Vertical lines indicate the value of F imp at the extrema 
(mg,g )-plane[(b)| 

V. RESULTS 



We calculate a perturbative series of the action (14.21) up to the fifth order. From this 
series we obtain the improved series according to the prescription presented in Section IHIl 
We identify plateaux as accumulations of extrema of observables with respect to artificial 
parameters. In the case of a model which exhibits phase transition, distribution of ex- 
trema and the value of observables in a plateau region change drastically around a phase 
transition point [^(J. Thus we examine the distributions of extrema of the improved series 
F imp (a, m 2 ; ao, m 2 ,) and A imp (a, m 2 ; ao, ttIq) with respect to ao and for each a and m 2 . 
In this paper we sweep a and m 2 from 0.0 to 10.0. This is because when these parameters 
are beyond 10.0 we do not observe a drastic change in the distributions of extrema and the 
value of observables around the extrema. 

In the following, we begin with describing detailed examinations at four characteristic 
points in the parameter space (a,m 2 ). 



(a } m 2 ) = (0.0,0.0) 

This corresponds to the pure Yang-Mills matrix model only with the four-point inter- 
action terms. 
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FIG. 2: Distribution of extrema of F imp (10.0, 0.0; olq, itiq) with respect to rriQ and Qo- Vertical lines 
indicate the value of F imp at the extrema (a) and its 2-d plot in (toq, ap)-plane |(b)| 

The distribution of extrema of F imp (ao, tuq) with respect to tuq, «o and its value at each 
extremum are shown in Fig. HJ We can find two accumulations of extrema as shown 
by the dashed circles A, B. These two accumulations also appear for A imp (ao, wig) and 
the improved series for the expectation value of the second moment of eigenvalues 

C( a V)4W^tr(4)\. (5.1) 



3 ^ \N 

n 

Except for these accumulations, there are gentle plateaux in the regions where the 
value of | mo | is very large. However, we do not take into account these flat regions 
as plateaux because these are nothing but asymptotic behaviors. These asymptotic 
behaviors always appear and we neglect them hereafter. By comparing the values of 
the improved free energy, we regard the plateau of region A as the true vacuum. The 
concrete values of the extrema and those of F irap (ao, m^) there are listed in Table [J 
and those for A imp (ao, ml) are listed in Table HT1 Because A imp is zero (Table IITI). this 
vacuum can be recognized as the Yang-Mills vacuum. 

Thus this parameter region corresponds to the Yang-Mills phase. 
(a,m 2 ) = (10.0,0.0) 

The distribution of extrema of F imp (a ,mQ) with respect to m , a and its value at 



each extremum are shown in Fig. [2j From the Fig. 2(b) we find an accumulation of 
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extrema. We point out that this accumulation is divided into two sets A and B by 
examining _F imp . As for the region A, the values of F imp at extrema go to large negative 
values as the order increases. On the other hand the values of F imp at extrema in the 
region B remain 0(1). They also appear for A imp . The concrete values of the extrema 
and those of F imp (a ,mg) there are listed in Table [TTT| and those for A 1Tap (aQ,mfy are 
listed in Table |IVl 

One might think that the region A cannot be regarded as a plateau because the 



values of F imp are not stable against the increase of the order (Fig. 2(a)). These 
values are all beyond 0(1) and in particular, the values of _F imp and y4 imp at the 
fifth order reach to about O(10 3 ) and (9(10 4 ), respectively, as shown in Table II III 
and Table IIVI In general the improved perturbation theory is an approximation by 
polynomials of finite order and therefore the value of 0(N 2 ) which goes infinite as 
iV — > oo is expected to be approximated as large amount of value. As we go to higher 
orders, the approximation would become better and would provide a larger value in 
such a case. In this sense, instability of an approximated value against an increase 
of order is not always problematic. Thus we identify the region A as a plateau. By 
comparing the values of F imp in the region A and B we regard the plateau of the region 
A as the true vacuum. 

Because the value of A imp in this plateau is O(10 4 ) (Table HVT) . this vacuum can be 
recognized as the fuzzy sphere vacuum and this parameter region corresponds to the 
fuzzy sphere phase. 

(a,m 2 ) = (4.0,1.5) 

The distribution of extrema of F imp (a ,mg) with respect to m , «o and its value at 
each extremum are shown in Fig. [3j We can find two accumulations of extrema. They 
also appear for A imp . By comparing the values of _F imp , the region B is expected as the 
true vacuum. However the values of C imp at extrema in this plateau become negative 
(Fig. H]). Thus we determine this plateau is unphysical. Then the plateau in the 
region A is recognized as the true vacuum which has negative 0(1) values of F imp and 
A imp . The concrete values of the extrema and those of F imp (a 0; m o), C imp { a Oi m o)' 
and A imp (a ,m 2 ) ) are listed in Table PVT IVTl and IVIH respectively. 

Because the value of y4 imp in the region A is negative 0(1) value (Table IVIII) this 
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FIG. 3: Distribution of extrema of i ?imp (4.0, 1.5; an> m o) w hh respect to mg and an- Vertical lines 
indicate the value of F imp at the extrema (a) and its 2-d plot in (mj, ao)-plane |(b)| 
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FIG. 4: Distribution of extrema of C imp (4.0, 1.5; ao, tuq) with respect to rriQ and ao- Vertical lines 
indicate the value of C imp at the extrema (a) and its 2-d plot in (m^, ao)-plane (b) 



parameter region corresponds to the 5-brane phase. 



(a,m 2 ) = (2.0,10.0) 



The distribution of extrema of -F imp (a , m^) with respect to m , an and its value at 
each extremum are shown in Fig. We can find an accumulation of extrema. The 
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FIG. 5: Distribution of extrema of i ?imp (2.0, 10.0; olq, ttiq) with respect to rriQ and Qo- Vertical lines 
indicate the value of F imp at the extrema (a) and its 2-d plot in (toq, ap)-plane |(b)| 

concrete values of the extrema and those of F imp (a , wiq) there are listed in Table IVIlTl 
and those for A imp (a , ml) are listed in Table HX1 

Because the value of A imp in this plateau is around zero (Table IIXI) . this vacuum 
corresponds to the Yang-Mills vacuum. Thus this parameter region corresponds to 
the Yang-Mills phase. 

Now we proceed to discussion in the case of generic a and m 2 . As an illustration we plot 
the values of F imp and A imp at their extrema with positive 6 m\ as functions of a when m 2 
is zero (Fig. [6] and Fig. [7]). From these figures we can divide the set of extrema into two 
branches: one branch consists of extrema for which the values of F imp and A imp remain of 
(9(1) against the increase of a, while the other branch is composed of extrema for which the 
values of F imp and A imp decrease far more than of O(10) as a increases. This property can 



be seen clearly in Fig. 7(b) We call the former as branch B and the latter as branch A. In 
fact each branch corresponds to an accumulation of extrema in (cto, m^-plane. The branch 
A, B corresponds to the plateau of the region A, B respectively at the above mentioned case 
(a,m 2 ) = (10.0,0.0). These two branches join together for small a, and thus we find only 



6 More precisely, we chose the region TOq = ~ 20, «o = — 20 ~ 20. This is because there is no accumulation 
except for foregoing asymptotic behavior in the region of |ao|, |wq| > 20. 
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FIG. 6: Values of F imp (a, 0.0; ao, tUq) at extrema as functions of a. 
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FIG. 7: Values of ^4 imp (a, 0.0; ao, tuq) at extrema (a) and its absolute value in logarithmic scale 
|(b)| as functions of a. 

one accumulation at (a,m 2 ) = (0.0,0.0). 

By comparing the values of _F imp (Fig. [6]), we identify the plateau corresponding to the 
branch A as the true vacuum as in the case of (a, m 2 ) = (10.0, 0.0). We take the expectation 
value of the Chern-Simons term as an order parameter. By examining the value of y4 imp in 
this branch we can determine what phase appears at each (a,m 2 ). If it takes a value of 
0(N 2 ), of Oil), and around zero, it corresponds to the fuzzy sphere phase, 5-brane phase, 



and Yang-Mills phase, respectively. Thus from Fig. 7(b) we conclude that there is the Yang- 
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FIG. 8: Values of A imp (5.0, m 2 ; oiq, Wq) at extrema as functions of m 2 . 

Mills phase in the region < a < 2, the 5-brane phase in the region 2 < a < 4 and the 
fuzzy sphere phase in the region a > 4. We also note that if the set of extrema is divided 
into branches, it indicates that the behavior of the improved series and the distribution 
of extrema change drastically there. Then this would imply a phase transition there as 



mentioned in the beginning of this section. In fact, Fig. 7(b) shows this is the case with 
our model, namely around the phase transition points we can observe the branching of the 
extrema. 

In the same manner, we observe the behavior of the value at extrema of A imp along an 
arbitrary path in (a, m 2 )-plane. As an example, we show values of A imp (a, m 2 ; a > a o) & t 
extrema as functions of m 2 when a = 5.0 in Fig. [HJ From these investigations we can obtain 
a phase diagram of the massive YMCS matrix model in (a,m 2 )-plane (Fig [9]). 



VI. CONCLUSIONS AND DISCUSSIONS 



In this paper we applied the improved perturbation theory to the massive YMCS matrix 
model and calculated the free energy, and the expectation values of the second moment of 
eigenvalues and the Chern-Simons term up to the fifth order. Though the original pertur- 
bative series was obtained about the perturbative vacuum corresponding to the Yang-Mills 
vacuum, the improved perturbation theory reveals the information on the different vacua, 
namely the 5-brane vacuum and the fuzzy sphere vacuum. 

As a result, we obtained the phase diagram of this model. The results are very simi- 
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FIG. 9: A proposed phase diagram of the massive Yang-Mills Chern-Simons matrix model. 



lar to the previous work based on the one-loop calculation around the fuzzy sphere type 
background. Since the improved perturbation theory is expected to include nonperturbative 
effects, we can expect that a gauge symmetry is generated as a subgroup of the original 
U(N) group even nonperturbatively in the 5-brane phase and the fuzzy sphere phase. In 
particular, dynamical generation of a non-Abelian gauge group occurs nonperturbatively in 
the 5-brane phase. However the precise positions of the critical lines in our phase diagram 
are somewhat ambiguous. For example, in our phase diagram there exists the 5-brane phase 
for 2 < a < 4 when m 2 = 0, which apparently contradicts with the Monte Carlo result in 
Ref. [7[. This is because our improved perturbative series is constructed from the standard 
perturbative series of lower order. It is obvious that in order to observe the phase transition 
clearly, a lower order analysis is not sufficient. As we proceed to higher orders, this ambigu- 
ity would decrease and we would find the phase transitions between the three phases more 

n 

precisely. Then a precise comparison to the phase diagram proposed in Ref. p] makes sense 
and we can discuss to what extent the perturbative analysis is valid in various parameter 
regions. 

Another result of this paper is a confirmation of validity of the improved perturbation 



theory. It is accomplis 
Monte Carlo method 



led by comparing our result in Fig. [7] with the precedent result of the 



221 ] . Although we can not see clearly the first-order phase transition 
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between the Yang-Mills phase and the fuzzy sphere phase, our result indicates a similar 
parameter region (a cr iu C ai ~ C(l)) as a critical point to that of the Monte Carlo method. 

It is worthwhile to mention the accumulation which was regarded as an unphysical plateau 
when (a,m 2 ) = (4.0, 1.5). This accumulation of the extrema shows the lower value of the 
free energy than that of the other plateau. Therefore it would be the true vacuum if it 
corresponds to a physical state. It is interesting that it has a positive expectation value of 
the Chern-Simons term which is an opposite sign compared to that in the usual 5-brane or 
fuzzy sphere vacuum. We considered this plateau unphysical because the expectation value 
of the second moment of eigenvalues is negative there. However if it becomes positive when 
we go to a higher order or relax the S'0(3)-symmetric assumption (14.41) . then we would have 
a nontrivial phase and the phase diagram would change. We leave it to a future work. 
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order 


a 


m o 


^? im p 


3 


-4.68729 x 10 y 


-8.30606 x 10 Y 


26.1027 




4.68729 x 10 9 


-8.30606 x 10 7 


26.1027 




0.0 


-2.40953 


1.5702 




0.0 


2.40953 


1.5702 




-2.31723 


4.29037 


1.70658 




2.31723 


4.29037 


1.70658 


4 


0.341271 


2.59243 


1.54085 




-0.341271 


2.59243 


1.54085 




3.40785 


5.75606 


1.71942 




-3.40785 


5.75606 


1.71942 


5 


-2.46938 


-4.30854 


1.54941 




2.46938 


-4.30854 


1.54941 




-1.49032 


-3.65313 


1.55424 




1.49032 


-3.65313 


1.55424 




0.0 


-2.54153 


1.51756 




0.0 


2.54153 


1.51756 




-0.524011 


2.9491 


1.5305 




0.524011 


2.9491 


1.5305 




-4.51448 


7.27529 


1.73146 




4.51448 


7.27529 


1.73146 



TABLE I: The numerical data of extrema for F im P(0.0, 0.0; a , ml) 
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order 


a 


2 

m 




3 


0.0 


-3.08221 


0.0 




-0.138803 


-2.94881 


0.000320748 




0.138803 


-2.94881 


-0.000320748 




0.0 


-2.82843 


0.0 




0.0 


2.82843 


0.0 




0.0 


3.08221 


0.0 


4 


5.45504 


-6.42054 


0.116302 




-5.45504 


-6.42054 


-0.116302 




0.0 


-3.7272 


0.0 




-0.290656 


-3.458 


0.000734582 




0.290656 


-3.458 


-0.000734582 




0.629447 


-3.26122 


0.000832601 




-0.629447 


-3.26122 


-0.000832601 




0.0 


-3.24237 


0.0 




0.270866 


-2.89579 


0.00388426 




-0.270866 


-2.89579 


-0.00388426 




0.0 


-2.66364 


0.0 




0.0 


2.66364 


0.0 




0.0 


3.24237 


0.0 




0.0 


3.7272 


0.0 


5 


0.0 


-4.42101 


0.0 




0.445076 


-4.03236 


-0.000992029 




-0.445076 


-4.03236 


0.000992029 




0.0 


-3.75696 


0.0 




0.965698 


-3.72636 


0.000171623 




-0.965698 


-3.72636 


-0.000171623 




-1.00261 


-3.69771 


-0.000168316 




1.00261 


-3.69771 


0.000168316 




0.0787589 


2.77319 


-0.000243519 




-0.0787589 


2.77319 


0.000243519 




0.0 


3.75696 


0.0 




0.0 


4.42101 


0.0 



TABLE II: The numerical data of extrema for ,4 imp (0.0, 0.0; a , m§) 
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order 


a 


2 


^? im p 


3 


9.43649 


-7.60688 


3.07945 




9.04356 


-6.42828 


3.06715 




-4.43165 


-5.84631 


1.67693 




1.43896 


0.862241 


-49.4706 




0.876274 


1.6869 


-9.32175 




-10.2822 


9.45705 


2.64061 


4 


10.5828 


-9.79928 


3.08573 




-0.00542322 


-2.96341 


7.73904 




1.53846 


0.937571 


-243.211 




1.32581 


1.08363 


-169.926 




0.00542358 


2.95278 


-4.62466 




-16.3116 


13.527 


2.75701 


5 


11.5522 


-11.9498 


3.08677 




-9.84509 


-9.44094 


1.5811 




0.791191 


-2.37935 


-14.4279 




1.42496 


0.986792 


-762.155 




1.58502 


1.01745 


-835.662 




1.06366 


1.73127 


-116.83 




-1.32478 


4.47375 


-2.3464 




-22.4794 


17.6733 


2.82902 



TABLE III: The numerical data of extrema for F imp (10.0, 0.0; a , mg) 



order 


a 


m o 




3 


9.94959 


-7.34049 


0.0901421 




-0.00271556 


-2.95803 


1.23634 




1.40096 


0.995306 


-106.684 




0.00271758 


2.95805 


-1.23634 


4 


-22.1864 


-15.5176 


-0.0340077 




11.4184 


-9.11834 


0.0886916 




0.798254 


-2.2225 


-8.76388 




1.501 


0.998731 


-519.979 




1.12189 


1.54166 


-73.7818 




-1.65882 


4.95124 


-0.637283 


5 


12.4504 


-10.8361 


0.0890299 




8.65517 


-6.44966 


0.101811 




-1.90908 


-4.94843 


1.46692 




1.1106 


-2.05508 


40.6892 




1.51445 


1.01601 


-1385.51 




1.42806 


1.18735 


-1840.99 




0.703556 


2.33845 


-37.6504 




-3.67883 


7.04633 


-0.431045 



TABLE IV: The numerical data of extrema for ^ im P(10.0, 0.0; a , ml) 
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order 


a 


2 


pimp 


3 


9.61196 x 10 1U 


-6.49178 x 10 Y 


25.8232 




-5.7603 x 10 9 


-4.47955 x 10 6 


22.5561 




-1.92635 


-3.3072 


-0.277788 




2.01317 


-3.09678 


1.21161 




1.71943 


-2.35192 


1.1326 




2.17535 


2.25552 


1.94668 




0.966819 


2.85037 


2.20386 




-7.70211 


8.53538 


2.99592 




1.82382 x 10 12 


1.63535 x 10 8 


8.64457 


4 


2.66644 


-4.5446 


1.27361 




-0.00980172 


-2.38145 


1.97284 




1.7947 


2.42488 


1.90002 




2.27369 


2.50761 


1.86237 




0.0104127 


3.78695 


2.27765 




-11.785 


11.6617 


3.04714 


5 


3.01798 


-5.86682 


1.29304 




-4.23407 


-5.18175 


-0.613976 




0.45086 


-2.18319 


-0.0677834 




2.00342 


2.38336 


1.82213 




2.2562 


2.7082 


1.82699 




1.20384 


2.99017 


1.95765 




-1.03239 


4.86228 


2.3355 




-15.9168 


14.8553 


3.08096 




TABLE V: The numerical data 


of extrema for F imp (4.0, 1.5; a , 


ml) 
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order 


ao 


<> 




3 


3.39189 


-6.74048 


-0.48585 




-3.47241 


-5.02842 


-1.38562 




3.80442 


-4.2263 


-0.560315 




2.58809 


-4.03884 


-0.497245 




3.1278 


-3.50164 


-0.581923 




-1.621 


-3.47036 


-1.18014 




1.0354 


-2.03806 


-0.841414 




2.12207 


2.47803 


0.786263 




2.04249 


2.50236 


0.873765 




2.07634 


2.85198 


0.859968 




0.919071 


3.14954 


0.586119 




1.21074 


3.2149 


0.766484 




-0.859116 


5.0785 


0.514263 




-5.67532 


8.38804 


0.292278 




-11.7333 


13.9904 


0.273057 




3.87597 x 10 12 


4.55414 x 10 8 


1.73042 x 10~ 8 




-4.05022 x 10 12 


3.90586 x 10 9 


0.0 


4 


3.07677 


-5.41897 


-0.488433 




2.0976 


-2.75252 


-0.626686 




-0.0734577 


-2.68203 


-0.00223366 




1.25144 


-2.16141 


-0.5948 




2.12631 


2.58726 


0.829878 




1.79011 


2.68894 


0.811098 




0.0567346 


4.05283 


0.54443 




-8.67938 


11.1616 


0.280595 




-2.06462 x 10 11 


5.2415 x 10 7 


1.0194 x 10~ 6 




2.78288 x 10 11 


7.59018 x 10 7 


2.4174 x 10~ 7 


5 


3.39189 


-6.74048 


-0.48585 




-3.47241 


-5.02842 


-1.38562 




3.80442 


-4.2263 


-0.560315 




3.1278 


-3.50164 


-0.581923 




0.330493 


-2.35159 


-2.49497 




2.04249 


2.50236 


0.873765 




2.07634 


2.85198 


0.859968 




1.21074 


3.2149 


0.766484 




-0.859116 


5.0785 


0.514263 




-11.7333 


13.9904 


0.273057 




-1.2577 x 10 11 


4.10542 x 10 7 


2.22914 x 10~ 6 




5.85754 x 10 11 


1.06032 x 10 8 


1.83291 x 10" 6 



TABLE VI: The numerical data of extrema for C imp (4.0, 1.5; a , mg) 
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order 


"0 


2 


^imp 


3 


2.34863 


-2.98275 


0.344616 




-0.00869776 


-2.30156 


1.30642 




1.98521 


2.40399 


-0.679792 




0.00529296 


3.80166 


-0.198496 


4 


-8.43836 


-7.87598 


-0.122983 




3.58394 


-4.43987 


0.312971 




0.553466 


-2.13643 


-0.541593 




2.11319 


2.50874 


-0.812064 




1.33838 


2.85121 


-0.576061 




-1.28207 


5.19139 


-0.156768 


5 


4.11377 


-5.56741 


0.310493 




-1.1206 


-3.34363 


2.1189 




0.722978 


-2.04584 


2.72295 




1.04265 


-2.01922 


0.734287 




2.08655 


2.54758 


-0.866595 




1.93161 


2.68347 


-0.877446 




0.744671 


3.4634 


-0.504538 




-2.62651 


6.65831 


-0.133145 



TABLE VII: The numerical data of extrema for yl imp (4.0, 1.5; a , ml) 



order 


a 


m o 


pimp 


3 


1.64718 


10.4191 


4.97903 




17.4705 


17.5462 


4.99209 




-24.2305 


20.3698 


5.03619 


4 


1.91975 


10.5121 


4.97903 




0.00745054 


10.8081 


4.97907 




0.0074495 


10.8081 


4.97907 




26.7536 


22.5472 


5.00089 




-36.51 


26.4821 


5.05249 


5 


5.0158 


2.06827 


-59.5252 




-5.75065 


2.28982 


-50.2504 




-12.1958 


5.41985 


9.34926 




12.4862 


6.04493 


6.64019 




1.76262 


10.5456 


4.97903 




1.95766 


10.6191 


4.97903 




-2.77239 


12.0802 


4.97928 




35.9114 


27.6808 


5.00761 




-48.7488 


32.7442 


5.06396 



TABLE VIII: The numerical data of extrema for F imp (2.0, 10.0; a , mg) 
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order 


a 


2 




3 


0.566682 


-0.941477 


48.9325 




0.460771 


-0.890109 


48.8427 




-0.043734 


-0.807144 


68.646 




0.113023 


-0.728224 


58.1744 




1.92014 


10.6554 


-0.00337663 




0.00372349 


10.8095 


-0.00330928 


4 


-1.33655 


-1.39229 


-116.063 




-0.166586 


-1.07679 


202.618 




-0.458842 


-0.94533 


57.275 




0.131775 


-0.934978 


108.256 




-0.128625 


-0.778924 


-145.652 




0.257705 


-0.733339 


974.75 




0.12539 


-0.674379 


885.926 




1.83575 


10.5838 


-0.00337685 




1.90363 


10.8436 


-0.00337683 




-3.28575 


12.5122 


-0.00307347 


5 


1.08531 


-1.82014 


59.1655 




0.511521 


-1.44973 


41.398 




-0.271338 


-1.42093 


434.787 




0.174771 


-1.23817 


73.8638 




1.90812 


10.5967 


-0.00337686 




1.33656 


10.6194 


-0.0033768 




1.87274 


11.0713 


-0.00337686 




-6.75945 


14.7516 


-0.00278608 



TABLE IX: The numerical data of extrema for ^ im P(2.0, 10.0; a , m§) 
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